Dynamical mass generation by source inversion: 
Calculating the mass gap of the Gross- Neveu model. 

K. Van Acoleyen*and H. Verschelde 
February 1, 2008 



Department of Mathematical Physics and Astronomy, 
University of Ghent, 
Krijgslaan 281,9000 Ghent,Belgium 

Abstract 

We probe the U(N) Gross-Neveu model with a source-term J^^. We find 
an expression for the renormalization scheme and scale invariant source J, 
as a function of the generated mass gap. The expansion of this function is 
organized in such a way that all scheme and scale dependence is reduced to 
one single parameter d. We get a non-perturbative mass gap as the solution 
of J = 0. In one loop we find that any physical choice for d gives good results 
for high values of N. In two loops we can determine d self-consistently by 
the principle of minimal sensitivity and find remarkably accurate results for 
iV > 2. 
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1 Introduction 



Asymptotic free massless field theories suffer from I.R.-renormalons, which 
signal that one is expanding around a wrong vacuum. To cure the theory 
from these renormalons the particles must acquire a mass. The only di- 
mensionful parameter available is A-^g^. A problem arises: in conventional 
perturbation theory A-^ appears only as a log not as a power. 

In the recent past, two new solutions, besides the 1/A'^-expansion, have 
been formulated for this problem in the context of the Gross-Neveu model 
[^, which is a very appealing model to test non-perturbative methods be- 
cause the exact result for the mass gap [|| is known. One of the two solutions 
consisted of a renormalizable version of the Optimized Expansion Q, while 
the other, found by one of us, presented a fully renormalizable action for 
the local composite operator Both methods require some tricks 

to maintain renormalizability. 

In this paper we present a third method which has the advantage of a 
straightforward renormalization. The idea is very simple: we add a source 
term to the Gross-Neveu Lagrangian, then we calculate the mass gap 
to obtain the relation m{J). After inversion^, we find an expression of J, 
the renormalization scheme and scale independent quantity associated with 
J, as a function of m. Putting J(m) equal to zero, gives us all possible 
solutions for the mass gap. Furthermore we will organize the expansion of 
J(m) in such a way, that all the scheme and scale dependence reduces to 
one single parameter d. This parameter is fixed by the principle of minimal 
sensitivity. 

In section ^, the method is explained in detail. Section ^, contains the 
numerical results of the one and two-loop calculations. 

2 The Method 

Let us perturb the U{N) Gross-Neveu model in 2 — e dimensions with a 
composite operator. 

C = md^jf^^if + ^fi'g'^{^^f - J^^ + C^t. (1) 

We have N Dirac fermions which enjoy a manifest U (A^)-symmetry. This 
theory is asymptotically free and for J = 0, possesses a discrete 75 invari- 

^The idea of inversion is not new, it was already formulated by Fukuda[^. He consid- 
ered the inversion of $(J), where "1> is a composite operator. An extensive review of the 
method and its applications can be found in M. 
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ance which in perturbation theory leads to a vanishing mass gap. It is also 
renormalizable with the counterterms: 



C^X. = i5z^9^7A.^ + \fi'6z,g^{^^f - 6z,J^^. (2) 

These counterterms depend on the scale /i and on the renormalization scheme. 
Because we are interested in the limit J — > 0, we will restrict ourselves to 
the mass independent renormalization schemes. In the context of the mass 
gap, every (mass independent) R.S., can be obtained from another (the 
MS-scheme for instance) by the following transformations: 

7= J(l + ao5^ + ai/ + ...), (3) 
f = g\l + bo9^ + b,g^ + ...), (4) 

ai,hi are finite numbers and J = J{fi),J = J{fi),g'^ = g^{fi),g^ = g^{fj-) is 
understood. From now on we take the convention that J,g'^, Pi,^,. . . , refer 
to the MS scheme, while J, g'^ , refer to an arbitrary R.S.. The scale 
dependence of J and g^ is governed by the 7 and (3- functions: 

d 

f^-Qj^J\jo,go,e = -Ji{g^) = -J{iog^ + iig^ + i2g^ + ■■■), (5) 

l^^gX,e^Pig^) = -2(/3o/ + A/ + /525' + ...)• (6) 

The R.S.-dependence of 7i,72,/32 (70 = To, /3o = Po^Pi = Pi) can be easily 
calculated using (P)-(|6|). In the formulation of our method we were led by 
the requirement that physical results must be independent of the scale and 
the scheme. 

We now define the scale and scheme independent J: 

J^Jfig^), (7) 

with 

^^§|^fi9')\9o,e = l{g')f{g')■ (8) 

The scale independence is seen immediately. As a solution of (ph, we find: 



2n X -i^n 1/71 7o/3l^ 2 , / 72 , 7i/3i 7o ,Pi.2 , 7o^2 

7i 7i70/?i I I X 
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with A an integration constant. We choose A = 1, any other choice imphes 
a different J : J(A) = AJ(0). This is not a problem because we wih look for 
solutions of J = 0. One can now easily check the scheme independence of 
J. 

In order to get J(m), we will first have to calculate m( J). A perturbative 
calculation of the mass gap generated by (|l|) will give us: 



m{J) = J{l+g\X', + X[\n^-^)+g\Y^ + Yi\n^ + Yi{}n^-^f) + . 
After inversion, we arrive at a relation of the following form: 

,2 



(10) 



J(m) = m(l + g''{Xo + Xi In 



m 



+g\Y, + In ^ + Y^iln ^ f) + ...)• 



(11) 



The coefficients which multiply a log can easily be related to the other 
coefficients and the coefficients of the 7 and /5 functions. Demanding @ 
and (Ml, we obtain: 



J(m) = m(l + g\Xo + ^ In — ) + g\Yo + + - Po)Xo) In — 

t2 



i.i 

4^ 2 



+7(^-/5o7o)(ln^)2) + ...). (12) 



Plugging this result together with (|9|), in (^ we finally arrive at: 



J 



_ 20. _ 
m{g 1^0 ~ 



1+5'Ui 



+9 i^i^ 



PI 

7i7o/?i 7o^:^ 
4/33 8/3o 



/3o' 



1 .71 7o/9i x 70 . ^ 
2V0 /3o^^ 2'''^^ 

7l/?l Jt)_/^^2 , 7o/?2 

4/3oVo 



72 

4/3o 4/32 



7? 



1 , 7oPi , ^ , /7o. 



7o.7o/3i 71 N , 71 



70 



Pi 



/3o)Xo In 



+ 



1/7| 
4 V 2 



Polo ) (In 



m2 2 



+ ... 

,2 



(13) 



m{g--o) 1 + ^2(^0 + ^1 In ^) + <7'(5o + ^i In ^ 
+S,(ln^)2) + .. 



(14) 
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Due to the scheme independence of J, the coefficients Ai, Bi are independent 
of the mass renormahzation (|3|) . They depend only on the couphng constant 
renormahzation (^. We find ^ 

^0 = ^-&o^ 
Ai = M 

So = B-o + ]^{^ + l)b^-^b^ + boMl-^) (15) 

^Po 

So the expansion of J in powers of is still highly scheme {g'^,bo,bi) and 
scale (/u) dependent. We now show that this dependence is reduced to one 
simple parameter d, by using another expansion. Following an argument 
of Grunberg [^(II,A), we know that J, a scheme and scale- invariant quan- 
tity, depends only on m and Ajjg : J = mFij^). Reorganizing (|l^) as 

an expansion in powers of — 3— will leave one arbitrariness: we can 



/Join 



A=7 

MS 



equally well expand in ^ . To show this explicitly we will need 

MS 

some formulas: 

/3iln(ln^) 




2 
A2 



l3o 



In(ln^) -ln(ln^) + 



(/3oln 



/3oln 



^2 



(16) 



A = A^exph^] (see e.g. [9]) (17) 
^Po 

P2 = {bl-b,)Po + Pibo + P^ (seei). (18) 

Using these formulas^ together with (14) and (|l5|), we find the master for- 

2 The fastest way to arrive at (|l5|), is the substitution of (^) in formula ([l4|), written in 
the MS-scheme. Alternatively one can also use the scheme-dependence of 7i,72,/32,^o 
and Yq, checking explicitly the mass- renormahzation independence. 

^Note that the expansion ( [l^ is only valid for large positive values of /3oln 
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mula: 
J = 



m{Po\n- ^ + d)''i^o X 



A 2 



1 + 



(/3oln 



A 2 



+ d) 



2/^0 AlTTc PO 



2/?o 



(/3oln-^+d)2 



2P0 PO Ajg^ Po 



+ (|)X[M.n 



, m2 70 . 70 In , 70 , rn? 



'■MS 



MS 



2/3oVo /3o' 



.2^ 70 / 70 
4/3o^2/?o 

7o/3i , 7o/3i 



+ 



2/3o' 2/3o^ 

3 



Atto Po ^P( 



^M5 



(19) 



MS 

m? 



with d = 6o — /3o In 

We can now recover the original Gross-Neveu model, by putting Jq, the 
naked source, equal to zero: 



Jo{m) ~ J{m) = 
One finds a non-perturbative mass gap, that is a solution of 

m , , „ , ,s ^ r. 1 



(20) 



^A- 



MS 



m ^0 
(/3oln-— 2 +ci)^ 

^MS 



1 + 



(/3oln 



''■MS 



(/?oln^+d) 
+ .. 



= 



(21) 



The total series is of course cZ-independent, but one can only calculate it 
up to a certain order, this will give us a mass gap that depends on d. 
One can check that the d-dependence of the order n truncated series is 

O ^ ) . The only sensible thing one can now ask for d is 



MS 
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minimal sensitivity ^\d=do — 0; hope that the expansion parameter 
( ] is small enough to justify the truncation .RWe finally note 

MS 

that, as in one needs the (n + l)-loop divergencies (7,/?), to get an 

n-loop result. 

3 Numerical Results 

The exact expression for the mass gap was obtained in 0]: 



1 

r(i~A)^ 



^=i^e)''—. -^MS' (22) 



where A = 2N-2 • Expanding in powers of N we have: 

m = il + ^(i^^^f^) + 0{j,))Aj^. (23) 

Comparison of our results with the exact result will give us an explicit check 
of our method. 

The 2-loop result for m(J) in the MS-scheme was given in [^(section 

4): 

,n 2r^-V2i J\ 4.(iV-l/2)(iV-3/4) 2 

(iV-l/2)(jV-l/4) J2 jV-1/2 
+ ^2 + ^^(0.737775 - -)]J(24) 

After inversion one finds: 

2 2 4 1 2 

Jim) = m[l + |^[(iV-l/2)ln-^] + |^[-(iV-l/2)(ln-^)2 

--{N - 1/2) In - 4{N - l/2)(0.737775 - -)]] . (25) 



So we have: 



Xo = 

7r2 ^ 96 



Yo = Z2"^(0-737775 - — ). (26) 



^We cissume that the series is at least asymptotic and that the size of the terms is more 
or less determined by the size of the expansion parameter. 
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The (3 and 7-functions have been calculated (in the MS-scheme) up to three 

/9i = -^7z^> P2 = -- — (27) 



loops by Luperini and Rossi |1C] and Gracey and Bennett ||Tl|, [|l2|, 

N-l „ iV-1 _ {N -l)[N -7/2) 



70 



N -1/2 



vr 



71 



47r2 
TV - 1/2 



72 



167r3 

(iV- l/2)(Af -3/4) 

47^3 



■ (28) 



One can now easily check (12) 



3.1 one-loop results 

The 1-loop result for the mass gap is found as the solution of the 1-loop 
truncation of J = 0. Putting the numbers in (19), one finds the mass gap 
as the solution of the transcendental equation: 



1 



Af-l 



In- 



+ dn 



2N -I ^ , m 

l + 47r(i + 21n[ln 



8(A^- 1 



2d7r 



A 2 ^ 



1 



0(29) 



First, the bad news: there exists no do for which the minimum-sensitivity 
criterium holds. Indeed, taking the derivative to d of ( p9[ ) gives us an equa- 

0: 



tion which must also hold if ^ 



TT 



'Af-1 



In- 



+ dvr)^ 



2N 



8{N 



1 + 47rd+ 21n[ln 



m 



A 2 



2diT 
N- 



TT 



N-l 



In 



+ dn 



2N 



8{N 



N-l 



In 



+ d-K 



(30) 



or using (I 



N-l 



In 



+ dTT {2N - 1) 



(31) 



So we would find a negative expansion-parameter, which is clearly not con- 
sistent with (p^), since the last term gets an imaginary value. 

Now the good news: for any physical choice for d we find reasonable re- 
sults. Furthermore we can show that any physical d gives the exact N ^ oo 
limit. This happens because every natural coupling constant renormaliza- 
tion (0) leads to a value for bo ^=°° with a dependent on the condition, 
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N 


mil 


mi2 


i V LaJ 


1 /N 


2 


210.7% 


268.9% 


-46.3% 


-21.9% 


3 


31.9% 


40.6% 


-32.5% 


-12.2% 


4 


14.7% 


19.0% 


-24.2% 


-7.0% 


5 


9.1% 


11.9% 


-19.1% 


-4.5% 


6 


6.5% 


8.5% 


-15.8% 


-3.1% 


7 


5.0% 


6.6% 


-13.5% 


-2.3% 


8 


4.0% 


5.4% 


-11.7% 


-1.8% 


9 


3.4% 


4.5% 


-10.4% 


-1.4% 


10 


2.9% 


3.9% 


-9.3% 


-1.1% 



Table 1: one loop results 



defining the R.S.. So if we take = m in (|l^), we find d ^=°° and 
equation (^) becomes: 

N^oo 2a 

1 = 5 , or m = AiTc. 32 

ln-^ + 2a 

We now calculate the mass gap, using two different conditions for the 
coupling constant renormalization. The conditions are defined on the 1 
P.I. 4-point function. In zero loops this originates from the interaction- 
term ^{^^) . In one loop, one finds, after a straightforward calcula- 



2 2 

tion, an effective interaction of the form: ^(^'^') 1 + g'^f{J,s,t,u,fi) 



+ 



9\ 

ble natural conditions are: 



other operators 



, where s, t, u are the Mandelstam- variables. Two possi- 



f{J,s,t,u,^)\s^t=u=o = 0, (33) 
ii=J 

/(J,s,t,u,/x)L^,__,2 = 0. (34) 

(U=J 



They lead to two values for d, di = ^ ^^'^ and ^2 = — ^ ^ — V5ziL_ jf 

we now use these values in (p9|), we find two solutions mn and mi2- 

The deviation from the exact result for the N ^ oo limit, the first order 
in jj, and our two one- loop solutions mu and mi2 have been displayed in 



( 2 Vs'^'^Vs-i 



table 3.1 in terms of percentage. Our one-loop results are acceptable for 
> 4. They lie somewhere between the N = oo and 1/A^ approximations. 
For iV = 2, 3 the results are bad. We could try to find another renormal- 
ization condition, which produces better results, but this would still leave 



9 



M 


mi 




i/^Po in[m /A-jg^j -|- aojTT 


o 

I 


±Z070 


/ 
/ 


/ 

/ 


Q 
O 


n Q% 

u.y /o 


1 zL 




4 


-1.0% 


2.1 


0.34 


5 


-1.5% 


2.8 


0.28 


6 


-1.6% 


3.4 


0.24 


7 


-1.6% 


4.1 


0.20 


8 


-1.5% 


4.8 


0.18 


9 


-1.4% 


5.5 


0.16 


10 


-1.3% 


6.1 


0.14 


15 


-1.0% 


9.5 


0.10 


20 


-0.8% 


12.8 


0.07 



Table 2: two loop results 



us dissatisfied. After all, the value for d is determined by some external 
physical (what does that mean anyway?) condition and does not result in 
a self-consistent way from the calculations. Fortunately, this changes with 
the two-loop calculations. 

3.2 two-loop results 

Putting the numbers in (|l9|) and equating it to zero, we find, for > 2, 
one and only one value for d (and consequently for m), determined by the 
minimal sensitivity condition ^ = 0. In table |3.2| we display the deviation 
(in terms of percentage) from the exact result, the value for d and the value 
for the expansion parameter. We find excellent (< 2%) agreement. 

The N = 2 case is special, one finds no minimum. The graph for m(d), 
is displayed in figures [l| and ^. There is a gap with no solution for m in 
the region where the minimum is expected. With the minimal sensitivity 
prescription in mind, we can still estimate the mass in the range of ±20% 
deviation. The generic N = 3 case is plotted in figure I for comparison. 
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4 Conclusions 

In this paper we have presented a general method for dynamical mass gen- 
eration in asymptotically free theories. When tested on the Gross-Neveu 
model, remarkably accurate results were found for the two-loop calculation. 

Let us now try to recapture the essence of the method and its connection 
with the sign of /3o- The general idea was to probe the theory with a R.S.- 
and scale invariant source J. Without the inversion one would arrive at a 
relation of the following form: 



m( J) = J(/3o In 



A 2 



70 

+ d) 2/3o 



1 + 



1 



(/3oln. 



J2 



+ d) 



(/3oln-^ + d)2 



+ 



(35) 



The important point is that this expansion is only valid for large positive 



75 



values of (/jp In m . + d). For an I.R.-free theory (/3o < 0) this means 



< AjTo^ exp 



.2 3n 



'^0 , so we can take the limit J — s- and one finds no mass 
gap. For asymptotic free theories, the expansion ( |35j ) is invalid for small 
J, so it cannot be used to probe the theory around J = 0. The expansion 



of the inverse function J(m) (19) on the contrary, remains well defined in 
the limit J — > 0, provided that a positive solution for m exists and that the 

2 

corresponding expansionparameter l/jPo In ."^ ^ -|- d) is not too big. We 
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found that this was the case for the Gross-Neveu model if we fixed d by the 
minimal sensitivity prescription. Calculations on other models and on QCD 
are in progress. 
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